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ABSTRACT 

In this work, we consider the stochastic background of gravitational waves 
(SBGWs) produced by pre-galactic stars, which form black holes in scenarios of struc- 
ture formation. The calculation is performed in the framework of hierarchical structure 
formation using a Press-Schechter-like formalism. Our model reproduces the observed 
star formation rate at redshifts z < 6.5. The signal predicted in this work is below 
the sensitivity of the first generation of detectors but could be detectable by the next 
generation of ground-based interferometers. Specifically, correlating two coincident ad- 
vanced LIGO detectors (LIGO III interferometers), the expected signal-to-noise-ratio 
(S/N) could be as high as 90 (10) for stars forming at redshift z ~ 20 with a Salpeter 
initial mass function with slope x — 0.35 (1.35), and if the efficiency of generation of 
gravitational waves, namely, ecw is close to the maximum value ~ 7 x 10^"*. However, 
the sensitivity of the future third generation of detectors as, for example, the European 
antenna EGO could be high enough to produce (S/N) > 3 same with eqw ~ 2 x 10^^. 
We also discuss what astrophysical information could be derived from a positive (or 
even negative) detection of the SBGWs investigated here. 
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1 INTRODUCTION 

Gravitational waves (GWs) are a natural consequence of 
Einstein's theory of general relativity (GR). GWs will open 
a new astronomical window for tiie study of the Uni- 
verse transforming the research in GR into an observa- 
tional/theoretical study. In particular, the opening of the 
full electromagnetic spectrum to astronomical observation 
during the last century expanded our comprehension of the 
Universe. In this century, observations across the gravita- 
tional wave spectrum will provide a wealth of new knowl- 
edge, including the possibility of studying the period when 
the first stars were formed in the Universe in the end of the 
so-called 'dark ages'. 

The information provided by GWs is different when 
compared to that provided by electromagnetic waves. GWs 
carry detailed information on the coherent bulk motions of 
matter, such as those produced by the collapse of stellar 
cores generating, for example, black hole remnants. On the 
other hand, electromagnetic waves are usually an incoherent 
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superposition of emissions from individual atoms, molecules, 
and charged particles. 

Because of the fact that GWs are produced by a large 
variety of astrophysical sources and cosmological phenom- 
ena, it is quite probable that the Universe is pervaded by a 
background of such waves. Collapse of Population II and III 
stars, phase transitions in the early Universe, cosmic strings, 
and a variety of binary stars are examples of sources that 
could produce such a putative background of GWs (see, e.g., 
Maggiore||2000[ |de Araujo, Miranda fc Aguiar||2002l [20041 



ISandick et al.|2006||Suwa et al.|2007||Giovannini|2009| among 
others) . 

Note that the indirect evidence for the existence of 
gravitational waves came first from observations of the or- 
bital decay of the Hulse-Taylor binary pulsar (IHulse & Tay- 



lor||1974 1975a|b I. Direct detection though and analysis of 



gravitational-wave sources are expected to provide a unique 
insight to one of the least understood of the fundamental 
forces ( Belczynski, Kalogera fc Bulik||2002 1. They will also 
allow us to investigate the physical properties of objects that 
do not emit any electromagnetic radiation as for example 
isolated black holes. 

A number of interferometers designed for gravitational 
wave detection are currently in operation, being developed. 



2 E. S. Pereira and 0. D. Miranda 



or planned. In particular, the high frequency part of the 
gravitational wave spectrum (lOHz < / < 10* Hz) IS open 
today through the pioneering efforts of the first-generation 
ground-based interferometers such as LIGO. While detec- 
tions from this first generation of detectors are likely to be 
rare, the advanced LIGO upgrade may detect, among oth- 
ers, the stochastic signal generated by a population of pre- 
galactic stars. 

Thus, in the future, it may be possible to use GWs 
as a tool for studying the star formation at high redshifts. 
In particular, from the theoretical point of view, it can be 
found in the literature several works discussing this possibil- 
ity. For example, the gravitational wave background (GWB) 
generated from the core collapse supernovae resulting in 
black holes at high redshifts has been discussed by |Fer 



rari, Matarrese fc Schneider] ( |1999| ); |de Araujo, Miranda fc 



Aguiar (2000 20041 among others. On the other hand, the 



calculation made specifically for Population III supernovae 
resulting in black holes is presented in |de Araujo, Miranda| 

fc Aguiar] ( [2OO2I ). 

More recently, Sandick et al. ( 2006 1 calculated the 
GWB from Population III stars with the cosmic star for- 
mation history in the framework of hierarchical structure 



formation. On the other hand, Suwa et al. ( 2007 1 presented 



the GWB spectrum of Population III stars by calculating 
the GW waveforms based on results of hydrodynamic core- 
collapsed simulations (see also Suwa et al.|2007 l. It is worth 
stressing that in all of these works, one of the most impor- 
tant parameters responsible to characterize the GWB is the 
cosmic star formation rate (CSFR). 

Goncerning to the CSFR at high redshift, our knowledge 
is mainly based on numerical simulations performed by hy- 
drodynamical codes in a A-CDM cosmology. Certainly, these 
simulations must reproduce the observable Universe at red- 
shifts z < 6. In particular, the evidence for the existence of 
a large star formation at high redshift comes from, among 



others, the Gunn-Peterson effect ( Gunn & Peterson 1965 1 

~^Zp, found in high — 2; Lya 



and from the metallicity of ~ 10 z/© 
forest clouds ( [Songaila fc Gowie|1996{|Elhson et al.|2000 |. 

These results are consistent with a stellar population 
formed at z > 5 { Venkatesan|2000 1 . However, measuring the 
CSFR from observations requires a number of assumptions, 
with the form of the dust obscuration corrections and the 



stellar initial mass function (|Kroupa 2007 
tham, fc Hopkins|2008| 



Wilkins, Tren- 



Our main goal in the present paper is to discuss how 
the detection of a GWB could be used to give us some in- 
sight on the CSFR. This kind of study could also be used to 
constrain the fraction of massive stars that generates black 
holes at high redshift, and the efficiency of production of 
GWs by black holes whose distribution function is presently 
unknown. To do so, we use a hierarchical structure formation 
model similar to that developed by Daigne et al.| ( [2b 06). 

However, in our model the CSFR is obtained in a self- 
consistent way. That means, we solve the equation governing 
the total gas density taking into account the baryon accre- 
tion rate, treated as a infall term, and the lifetime of the 
stars formed in the dark halos. 

The paper is organized as follows. In Section 2, we 
present the Press-Schechter-like (PS) formalism used to de- 
termine the comoving abundance of collapsed dark matter 
halos. In Section 3, we discuss how to obtain the CSFR from 



the hierarchical model. In Section 4, we present the formal- 
ism used to characterize the GWB. Section 5 presents our 
conclusions. 



2 HIERARCHICAL FORMATION SCENARIO 

Press and Schechter (hereafter PS) heuristically derived a 
mass function for bound virialized objects in 1974 (Press 



fc Schechter 19741. The basic idea of the PS approach is 



define halos as concentrations of mass that have already left 
the linear regime by crossing the threshold 5c for non-linear 
collapse. Given a power spectrum and a window function, 
it should then be relatively straightforward to calculate the 
halo mass function as a function of the mass and redshift. 
However, it is worth stressing that the exact definition of 
the mass function, e.g., integrated versus difi'erential form or 
count versus number density, varies widely in the literature. 
To characterize different fits, it can be introduced the scale 
differential mass function f{a, z) ( Jenkis et al.|200T l defined 
as a fraction of the total mass per lna~^ that belongs to 
halos. That is, 



dp/pB _ M dn{M,z) 
din (7-1 ^ pb{z) dln[(7-i(A/,2)] 



(1) 



Where n(M, z) is the number density of halos with mass M, 
/5b (z) is the background density at redshift z, and a{M,z) 
is the variance of the linear density field. As pointed out by 



Jenkis et al. (2001 1, this definition of the mass function has 



the advantage that it does not explicitly depend on redshift, 
power spectrum, or cosmology; all of these are contained in 
a{M,z) (see also |Lukic et al.|2007 l. 

To calculate cr{M,z), the power spectrum P{k) is 
smoothed with a spherical top-hat filter function of ra- 
dius R, which on average encloses a mass M {R — 
[SM/ATTpBiz)]"-^^). Thus, 



a{M,z) = 



D^z] 
27r2 



k'^P{k)W^{k,M)dk, (2) 



where W{k, Al) is the top-hat filter: 



W{k, M) 



{kRf 



[sin(fci?) - kRcos(kR)] 



(3) 



and the redshift dependence enters only through the growth 
factor D{z). 
Then, 



a{M,z) = a{M,0)D{z). 



(4) 



In the more general case of a Universe with matter and 
a cosmological constant, the exact solution for the growth 



function is well approximated by (Carrol, Press & Turner 



19921 



D{a): 



5D,in(a) a, 



2[i ~ QAia) + nT + ^n^{a)] 



(5) 



where the relative density of the i— component is given by 
Qi — pi/pc and 'i' applying for baryons (b), dark energy (A), 
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and total matter (m), while a = 1/(1 + z) is the cosmological 
scale factor. 

As usual, the primordial power spectrum is assumed to 
have a power law dependence on scale, that is, P{k) oc k". 
For a scale- invariant spectrum the spectral index is n = 1. 
The rate at which fluctuations grow on different scales is 
determined by an interplay between self-gravitation, pres- 
sure support and damping processes. These effects lead to a 
modification of the form of the primordial power spectrum 
that is expressed in terms of a transfer function T{k). Thus, 
we have: 



dt 

dz 



9.78/i-iGyr 



(f + z)^Q,K + ^n.{l + Zf 



(12) 



In Eq. ( 10 1 we have used as upper limit Afn 



10 Mq . This choice permits a reasonable computational 
time to run the models. Moreover, models with Mmax = 
lO'^* Mq showed no considerable difference in the results. In 
the next Section, we discuss how to obtain the CSFR from 
the hierarchical scenario here described. 



P{k) = BkT{k), 



(6) 



where the normalization factor B is determined observation- 
ally. 

For the transfer function, we consider ( [Efstathiou, Bond 
fc White| 19921 ) : 



r(fc) 



(7) 



{1 + [ak + {hkf/^ + (ckYYy/" ' 

with u = 1.13, a = (6.4/r)/i"^Mpc, b = (3.0/r)/i"^Mpc, 
c = (1.7/r)/i-iMpc, and F = 9.^h e-"''''+^/"°') is the 



so-called shape parameter of the power spectrum (Bardeen 
et al.|1986||Peacock|199"9l ). 



We use throughout this work the mass function fit pro- 



posed by Sheth fc Tormen (19991. That is. 



/sT(cr) = 0.3222W — - exp 

V TT O" 



2o-2 



1 + 



a5? 



(8) 



where a — 0.707 and p — 0.3. 

At redshift z, the comoving density of dark matter halos 
in the mass range [M, M + dAI] is fsT{o')dAI , with (see, in 
particular, Daigne et al.|2006 l 



PDM = / fsT{cT)MdM, 
10 



(9) 



where pdm is the comoving dark matter density. 

We consider that the baryon distribution traces the 
dark matter distribution without bias. Thus, the density of 
baryons is proportional to the density of dark matter. The 
fact that stars can form only in structures that are suitably 
dense can be parameterized by the threshold mass Mmin. 
Thus, the fraction of baryons at redshift z that are in struc- 
tures is given by 



J A 



/sT (CT)MdM 



f ST (a) MdM 



(10) 



With this definition, the baryon accretion rate ah{t) 
which accounts for the increase in the fraction of baryons 
in structures is given by ( Daigne et al.||2006 l: 



ah{t) = r^bpc 



d~z) 



dMz) 



dz 



(11) 



where pc ~ SHq /8ttG is the critical density of the Universe. 

The age of the Universe that appears in ( |11[ ) is related 
to the redshift by: 



3 THE COSMIC STAR FORMATION 

In hierarchical models for galaxy formation the first star- 
forming halos are predicted to collapse at redshift z > 
20, having masses ~ IO'^Mq ( [Salvadori, Schneider fc Fer- 
2007j). In particular, the star formation history for a 



'galactic-like system' is determined by the interplay between 
incorporation of baryons into collapsed objetcs (stars, stellar 
remnants, and smaller objects) and return of baryons into 
diffuse state (gaseous clouds and intercloud medium of the 
system) . 

The later process can be two-fold: (a) mass return from 
stars to the 'interstellar medium of the system' through, for 
example, stellar winds, and supernovae, which happens at 
the local level; and (b) net global infall of baryons from out- 
side of the system. The former process is a well-known and 
firmly established part of the standard stellar evolution lore 
(see, e.g., [Chiosi fc Maeder||1986[ ), and although details of 
mass-loss in a particular stellar type may still be controver- 
sial, there is nothing controversial in the basic physics of this 
process. 

Thus, we use throughout this paper the basic process 
above described. To do that, we consider the baryon accre- 
tion rate ah{t), described by Eq. (Ill, as an infall term that 
supplies the reservoir represented by the halos. Therefore, 
the number of stars formed by unity of volume, mass and 
time is given by: 



d-^N 
dVdmdt 



= $(m)*(f), 



(13) 



where <l?(m) is the initial mass function (IMF) which gives 
the distribution function of stellar masses, and »I'(t) is the 
star formation rate. See that ^{t) is assumed to be indepen- 
dent of mass while "I>(m) is assumed to be independent of 
time. 

We use a Schmidt law ( |Schmidt|[T959l peS] ) for *(t). 
Therefore, 



d^M^, 
dVdt 



= *(t) = k[p^it)r 



(14) 



where is a constant that will be identified later, pg is the 



local gas density, and a — I. See that ( 141 shows that stars 



are formed by the gas contained in the halos. 

On the other hand, we assume that the IMF follows the 



Salpeter ( 1955 1 form 



$(m) — Am' 



-(l + x) 



(15) 
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where x — 1.35 (our fiducial value) and A is a normalization 
factor. 

The constant A is determined by the condition which 
all stars are formed into the mass range [minf,msup]. That 
is, 

ylm~<'+"'mdm = 1, (16) 

f 

and we consider minf — O.IM0 and msup — I4OM0 as limits 



in (16 1 



The mass ejected from stars, for example through winds 
and supernovae, is returned to the 'interstellar medium of 
the system'. Thus, we have: 



dVdt 



(m - mr)^{m)'^{t - Tm)dm, (17) 



m(t) 



where the lower limit of the integral, m{t), corresponds to 
the stellar mass whose lifetime is equal to t. In the inte- 
grand, rrir is the mass of the remnant, which depends on the 
progenitor mass, and the star formation rate is taken at the 
retarded time {t — Tm), where Tm is the lifetime of a star of 
mass m. 

For all stars formed in the halos, we use the metallicity- 



independent fit of Scalo (19861; Copi (19971 



logio(^n.) = 10.0-3.6 l0| 



M 



+ 



log 



M 

1° I m;; 



(18) 

where Tm is the stellar lifetime given in years. 

The mass of the remnant, rrir, in Eq. ( |17[ ) is calculated 
using the following assumptions: 

a) Stars with m < 1 M© have a high lifetime so they do not 
contribute for A/cj; 

b) Stars with 1 M© 5C m ^ 8 M© after evolving off the 
main sequence left carbon-oxygen white dwarfs as remnants, 
where 



: 0.1156 m + 0.4551; 



(19) 



c) Stars in the range 8 Mq < m 5C 10 M© after evolv- 
ing off the main sequence left oxygen-neon-magnesium white 
dwarfs with = 1.35 Mq; 

d) Stars with 10 Mq < m < 25 Mq explode as supernovae 
leaving neutron stars as remnants {rrir = 1.4 Mq); 

e) Stars with 25 Mq m 140 Mq produce black hole 
remnants. In this case, we consider that mr = Jnnc. Note 
that rriHe is the mass of the helium core before collapse (see 
|Heger fc Wo"osiey|[2002l ). Thus, 



13, 
24 ^ 



20 Up 



(20) 



Then, using equations (141 and (17 1 we can write an 



equation governing the total gas density (pg) in the halos. 
Namely, 



- dwt + dyi^+-(^)' (21) 

where ab(t), Eq. ( |11[ ), gives the rate at which the halos ac- 
crete mass. 



Numerical integration of (211 produces the function 
/9g(t) at each time t (or redshift 2). Once obtained Pg{t), 



we return to Eq. (141 in order to obtain the 'Cosmic Star 



Formation Rate' ^'(t). Just replacing ^'(t) by p* we can write 



(22) 



where the constant k represents the inverse of the timescale 
for star formation. Namely, k — 1/ts. 

We normalize the CSFR in order to produce — 
0.016 MQyr~^Mpc"^ at z = 0. With this normalization, 
we obtain a good agreement with both the present value 



of the CSFR derived by Springel fc Hernquist (20031, who 



employed hydrodynamic simulations of structure formation, 
and the observational points taken from [Hopkins] ( |2004[ 
2007|. 



The cosmological parameters used in our models are: 
Qa = 0.76, = 0.24, r2b = 0.04, o-g = 0.84, and Hubble 
constant Ho = 100/ikms"^ Mpc"^ with h = 0.73. 

In the Figure 1 we present the CSFR derived from our 
models in function of the threshold mass Mmin (see Eq. 10 1. 



We use a IMF with slope x = 1.35, Ts = 2.0 Gyr as timescale 
for star formation, and we consider that stars start to form at 
redshift Zini = 20. As can be seen, models with Mmin = 10® — 
IO^Mq have an excellent agreement with the observational 
CSFR at redshifts z < 6.5. See that the threshold mass 
Afmin act on the amplitude and the redshift (z*) at which 
the amplitude of the CSFR is maximum. 

The model with Mmin = 10^°Mq has a good agreement 
with data at z < 5. On the other hand, at more higher red- 
shifts (5 < 2 < 6.5) this model does not agree very well with 
the observational points. In Figure 1 we also included the 
CSFR derived by Springel & Hernquist (SH) for compari- 
son. Although our models with Mmin = 10® — IO^Mq have an 
amplitude greater than that derived by SH we can observe 
that both, SH and our models, fit very well the observational 
data. 

In Figure 2 we show the influence of Ts on the CSFR. 
We consider x — 1.35, Zini ~ 20, and we take A/min = IO^Mq 
for the threshold mass. Note that, Ts ^ 2.0 Gyr produces a 
gas comsuption timescale compatible with early type galax- 



ies ( de Freitas Pacheco|1997 K Thus, the first basic effect of 
increasing Ts is to shift the peak of the CSFR to lower red- 
shifts. That means, the higher the Ts parameter, the lower 
is the readshift where appears the peak of p* . In particular, 
the peak of p* is shifted from redshift 3.3 if Ts — 4.0 Gyr to 
6.1 if Ts = 1.0 Gyr. 

The parameter Ts is also related to the amplitude of 
Pi, (see also Eq. 22 1. See that considering Afmin = IO^Mq 
then the models with Ts = 2.0 — 3.0 Gyr are those that 
present the best concordance with the observational data. 
It is worth stressing that both parameters, A/min and Ts, 
produce similar effects on the results. That is, they act on 
the amplitude of p-^ and on the value of z*. In Figure 2 is 
also included the CSFR derived by Springel & Hernquist for 
comparison. 

In Figure 3 we see the influence of Zini on the evolution 
of pi,. The models have similar evolution at z < 5. How- 
ever, at larger redshifts, the model with Zini = 40 produces 
a CSFR higher than that obtained from Zini — 20. In partic- 
ular, the peak of the CSFR occurs at redshift 4.6 (5.5) for 
the model with zini = 20 (zini — 40). 
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Figure 1. The CSFR obtained from our models compared to 
the observational points (HP) taken from |Hopkins| l |2004] |2007| l . 
We used a standard Salpeter IMF (x = 1.35), and Tg = 2.0 Gyr 
as timescale for star formation. In this plot, we can see the in- 
fluence of Mmin (the threshold mass for halo formation) on the 
CSFR. The solid line represents the Springel & Hernquist CSFR 
(SH), the dashed line corresponds to Mmin = IO^Mq, the short 
dashed line corresponds to Mmin = 10^ M©, and the dotted line 
represents Mmin = IO^^Mq. 
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Figure 2. The influence of the timescale for star formation (ts) 
on the results. The solid line represents the Springel & Hern- 
quist CSFR (SH), the dashed line corresponds to Tb = 1.0 Gyr, 
the short dashed line corresponds to Ts = 2.0 Gyr, the dot- 
ted line corresponds to Tg = 3.0 Gyr, and the dot-dashed line 
represents Tb = 4.0 Gyr. These models have a threshold mass 
Mmin = IO^Mq, and a IMF with slope x = 1.35. HP stands for 
the observational CSFR | |Hopkins|2004l |2007| l . 



It is worth stressing that the CSFR is inferred from 
observations of the light emitted by stars at various wave- 
lengths. These observable samples are flux-limited, and thus 
the intrinsic luminosity of the faintest objects in the sample 
changes with redshift. This incompleteness of the samples is 
corrected by using a functional (Schechter function) to the 
luminosity function obtained from the observations them- 
selves. 

An important parameter on the determination of the 



Figure 3. Models with Tg = 2.0 Gyr and Mmin = IO^Mq but 
considering two different values for the initial redshift. The solid 
line corresponds to the Springel & Hernquist CSFR (SH), the 
dashed line corresponds to zi^i = 20, and the short dashed line 
corresponds to zi^i = 40. HP stands for the observational CSFR 
| |Hopkins|2004l|2007| |. 



CSFR is the obscuration by dust that is well known to affect 
measurements of galaxy luminosty at ultraviolet (UV) and 
optical wavelengths. Correcting for this effect is not always 
straightforward. Thus, there are large uncertainties associ- 
ated to the determination of the CSFR as can be seen from 
Figures 1 — 3 (see, in particular, Hopkins 2004 de Araujo 



Ifc Mirandal2005l who discuss these uncertainties with more 
details) . 



4 THE STOCHASTIC BACKGROUND OF 
GRAVITATIONAL WAVES 

In this section we use the CSFR (pi,) obtained from the 
hierarchical model to determine the stochastic background 
of gravitational waves (SBGWs) generated by stars which 
collapse to black holes. Initially, we present a quick overview 
on the formalism used to characterize a SBGWs because this 
subject is discussed in previous works (see, for example, de 

de Araujo fc 



Arauj o, Miran da fc Aguiar||2000| [2002] [2004] 



Miranda 2005; 



Miranda, de Araujo & Aguiar 



20041. After 



this quick overview we display and compare the results of 
the models considered. 

Let us write the specific flux received in GWs at the 
present epoch as 



dv 

And? dvohs 



dV, 



where 



dL„ 



(23) 



(24) 



is the comoving specific luminosity density (given, e.g, in 
ergs~^ Hz~^ Mpc~'^), which obviously refers to the source 
frame. See that dV is the comoving volume element, and di^ 
is the luminosity distance. 

The above equations are valid to estimate a stochastic 
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background radiation received on Earth independent of its 
origin. In the present paper li, can be written as follows 



dEaw 
dv 



Pt{z)^{m)dm, 



(25) 



where dEcw/dv is the specific energy of the source. Note 
that in the above equation Pi,{z) is the CSFR, and $(m) is 
the IMF. 

Thus, the flux F„(i^obs) received on Earth reads 



I'ohs 



1 dEow dv 
4:TTd? dp dvoh'. 



p^,{z)^{m)dmdV. (26) 



In particular, one can write the differential rate of pro- 
duction of GWs, for the case of a background produced by 
an ensamble of black holes, as follows 



dRBH = pi,^^^{m)dmdz. 



Using Eq.(27l it follows that 

1 dEaw dp 



47rdT dp duoh 



dRsu- 



(27) 



(28) 



Note that in the above equation, what multiplies di?BH 
is nothing but the specific energy fiux per unity frequency 
(in, e.g., ergcm"^ Hz~^), i.e.. 



/i/(l^obs) 



dEc 



dp 



i-Kd? dp dl'obs 



(29) 



On the other hand, the specific energy flux per unit 
frequency for GWs is given by ( Carr|1980 1 



r f ^ _ TTC^ ,2 

Jiy(PobB) — "^"-BH- 



(30) 



Also, the spectral energy density, the flux of GWs, re- 
ceived on Earth, F,^, in erg cm~'^ Hz~^ can be written 
as 



From the above equations one obtains 



_ 1 

"■BG — 



/ibhC^-Rbh. 



(31) 



(32) 



See that /ibh is the dimensionaless amplitude produced 
by the collapse of a star to form a black hole. Its expression 



is obtained from Thorne (19871. Thus 



.e„.r.4x,0-».«(^)(^)- 



(33) 



where eow is the efficiency of generation of GW's, and rrir 
is the mass of the black hole formed. 

It is worth mentioning that Eq. ( |33[ ) refers to the black 
hole 'ringing', which has to do with the de-excitation of the 
black hole quasi-normal modes. 

The collapse of a star to black hole produces a signal 
with frequency i^obs given by 



fobs 



1.3 X lO^'Hz 



(^)(l + .)-\ (34) 



where the factor (1 -I- z)~^ takes into account the redshift 
effect on the emission frequency. That is, a signal emitted 
at frequency Pc at redshift z is observed at frequency i/obs = 

p,{i + zy\ 

As discussed in the previous section, we consider that 
black holes are formed from stars with 25 Mq ^ m ^ 
140 Mq . The mass of the remnant is taken to be the mass 



of the helium core before collapse (see Eq. 20 1 



Another relevant physical quantity associated with the 
SBGWs is the closure energy density per logarithmic fre- 
quency span, which is given by 



1 d pew 



Pc d log //obs ' 

The above equation can be re-written as 



a 



GW = 



_ 47r 2 ,2 



fob: 
,3 „ 



(35) 



(36) 



'0 



Thus, given a star formation history, consisting of a 
star formation rate per comoving volume (CSFR), Pi,{z), 
and an initial mass function (IMF), <^{m), the stochastic 
background of gravitational waves produced by pre-galactic 
black holes can be characterize. 

Finally, to assess the detectability of a GW signal, one 
must evaluate the signal-to-noise ratio (S/N), which for a 
pair of interferometers is given by (see, for example, |Chris 



tensen|199 2"Flanagan'1993" Allen|T997l|de Araujo, Miranda 
fc Aguiar|2002, , 2004,; Rcgimbau fc de Freitas Pacheco|2006[ ) 



(S/N)2 



ml 

507r4 



T 



dp- 



l\p)nl^{p) 



"si'\f)si;"{p) 



(37) 



where is the spectral noise density, T is the integra- 
tion time, and "/{p) is the overlap reduction function, which 
depends on the relative positions and orientations of the two 
interferometers. For the "/{f) function we refer the reader to 
Flanagan ( |1993) who was the first to calculate a closed form 
for the LIGO observatories. 

Using the formalism described above and in the previ- 
ous sections we study a total number of 72 models varying 
the following parameters: 

a) the threshold mass (Afmin) for structure formation, where 
we consider the values IO^Mq, IO^Mq, and IO^^Mq; 

b) the exponent (x) of the IMF, where we consider x — 1.35 
('Salpeter exponent'), x = 0.35 which yields a higher number 
of black hole remnants than Salpeter IMF, and x = 2.35 
which produces a lower number of black hole remnants than 
Salpeter exponent; 

c) the timescale for star formation (r^), where we consider 
the values 1.0 Gyr, 2.0 Gyr, 3.0 Gyr, and 4.0 Gyr; 

d) the initial redshift (^ini) where star formation begins to 
occur. We take the values 20 and 40. 

On this set of models we use two criteria for selecting the 
best ones. The first criterion is to have good agreement with 
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Figure 4. The CSFR for models with M^in = W^Mq and good 
agreement with observational data. The main characteristics of 
these models are described in Table 1. 



Figure 5. The CSFR for models with M^in = lO^M© and good 
agreement with observational data. The main characteristics of 
these models are described in Table 2. 



Table 1. The main results of the models with Mmin = IO^Mq of 
the Figure 4. The signal-to-noise, (S/N), is presented for a pair of 
LIGO III (advanced configuration) interferometers. (S/N) is com- 
puted for one year of observation and we consider a gravitational 
wave efficiency eGW„iax = 7 X 10^*. 



Model 


^ini 


X (IMF) 


Ts Gyr 




(S/N) 


Al 


20 


1.35 


2.0 


4.6 


7.4 


A2 


20 


1.35 


3.0 


3.8 


3.8 


A3 


20 


0.35 


1.0 


4.4 


93.5 


A4 


40 


1.35 


2.0 


5.6 


9.8 


A5 


40 


1.35 


3.0 


4.6 


4.8 


A6 


40 


0.35 


1.0 


5.3 


119.9 



observational star formation data at redshifts z < 6.5[^ The 
second criterion is to produce a signal-to-noise (S/N) > 3 for 
a pair of 'advanced' interferometers. We consider this choice 
of (S/N) as reasonable for an adequate characterization of 
the SBGWs. 

Figure 4 presents the models with Mmin = lO^M© 
which satisfy the above criteria. 

Table 1 shows the main results for the six models 
Al — A6 which are presented in Figure 4. The efficiency of 
generation of GWs is taken from ^Stark fc Piran ( 1986 1 who 
simulated the axisymmetric collapse of a rotating star to 



black hole. We use their maximum value, namely, ecWmax — 
7 X 10~*. We will discuss below the dependence of eow on 
the results. 

See that to calculate de signal-to-noise ratio we consider 
that the integration time in Eq. (37 1 is one year. In the fifth 
column of Table 1 we present the redshift (z*) where the 
CSFR reaches its maximum value; in the sixth column we 
present the signal-to-noise ratio (S/N). 



^ We performed analysis over the models with (S/N) > 3. In 
particular, we determine the reduced chi-square defined as Xr = 
X^/doi (where "dof means "degrees of freedom"). We consider 
that models with Xr ^ 1 have good agreement with observational 
data. 



Note that there is possibility of detecting the SBGWs 
here proposed if eow is around the maximum value. Observe 
that for Salpeter IMF {x = 1.35) we obtain a significant 
(S/N) if Te ~ 2.0 - 3.0 Gyr. 

On the other hand, the models with x = 0.35 pro- 
duce the highest values for the (S/N). This happens because 
X = 0.35 produces a higher number of massive stars than the 
Salpeter IMF. In this case, the CSFR that fit the observa- 
tional data are those with Ts < 1.0 Gyr. See that the models 
with Ta < 1.0 Gyr have a short timescale for star forma- 
tion. These values for the parameter Ts are consistent with 
a high-mass stellar population. 

However, if the IMF of pre-galactic stars is close to 
X = 2.35 then there is no hope of detecting the SBGWs 
we proposed here, even for ideal orientation and locations 
of the LIGO interferometers. In particular, all models with 
X — 2.35 have (S/N) < 0.1, same for those models producing 
with excellent agreement with Hopkins data. Thus, the 
first conclusion is that it would be possible the detection of 
a background of pre-galactic black holes if the IMF of these 
objects is a:: > 1.35 and if ecw^ax ~ 7 x 10"'*. 

In order to see the influence of Mmin on the value of the 
signal-to-noise ratio we present in Figure 5 the models with 
Afmin = IO^Mq and that satisfy our two criteria as above 
defined. 

Table 2 shows the main results for the six models 
Bl — B6 which are presented in Figure 5. The first effect 
of Mmin is to shift 2* (for example, compare models Al and 
B2). That is, a halo with mass W^Mq collapses earlier than 
a halo with mass lO^M©. Thus, the maximum of star forma- 
tion for models with 10*Mq will be shifted to low redshifts. 

The second effect is on the amplitude of p* as discussed 
in the previous section. As the quantity of black holes is oc 
then increasing the value of Mmin the number of black holes 
formed will decrease. As a consequence, models with 1O*M0 
present a lower (S/N) than those with Mmin = W^Mq. 

The third effect can be seen comparing Tables 1 and 2. 
The models which satisfy the selection criteria with Mmin ~ 
IO'^Mq are those with Ts ~ 2.0-3.0 Gyr for x = 1.35. Other- 
wise, with Mmin ~ IO^Mq the selection criteria are satisfied 
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Table 2. The main results of the models with Mjr 



IO^Mq. 



Model 


^ini 


X (IMF) 


Tb Gyr 


2* 


(S/N) 


Bl 


20 


1.35 


1.0 


5.1 


11.9 


B2 


20 


1.35 


2.0 


3.8 


5.7 


B3 


20 


0.35 


1.0 


3.6 


72.8 


B4 


40 


1.35 


1.0 


5.2 


13.2 


B5 


40 


1.35 


2.0 


3.9 


6.2 


B6 


40 


0.35 


1.0 


3.8 


77.6 



CO 
I 



o 

Q. 



© 



0.1 



0.01 



0.001 




Figure 6. The CSFR for models with M^in = IQIOMq and good 
agreement with observational data. The main characteristics of 
these models are described in Table 3. 



if Ts ^ 1.0 — 2.0 Gyr for x = 1.35. This result can be un- 
derstood remembering that Ts also acts on the amplitude of 

That means, if we decrease the value of Ts the amplitude 
of increases (see, for an instance. Figure 2 and Eq. 22 1. 
On the other hand, as above discussed, if we increase the 
parameter Mmin, the amplitude of is reduced. Thus, if we 



change Mmin from 10 Mq to 10 Mq, we have to decrease 
the parameter Ts in order to obtain with good agreement 
with the observational data and also to produce (S/N) > 3. 

In Figure 6 we present the models with A/min = 1O^''M0 
(see details of the models in Table 3). Only those with 
X = 1.35 and Ts — 1.0 Gyr have a good agreement with 
observational data and produce (S/N) > 3. See that the 
difference between models CI and C2 is very subtle. 

This happens because the fraction of baryons in strcu- 
tures with M > 10^°Mq is very smaU at redshifts 20 - 40. 
Thus, Zini does not have strong influence on the evolution of 
the models CI and C2 at low redshifts. 

Figure 7 shows the density parameter Qgw as a func- 



tion of the observed frequency fobs- The density parameter 
increases at low frequencies and it reaches a maximum am- 
plitude of about 9.0 X 10~^ around 200 Hz in the model A6. 
On the other hand, model A2 [preaches a maximum ampli- 
tude of 4.2 X 10'** also around 200 Hz. See that both the 
maximum amplitude of Qqw and the high-frequency part of 
the spectra |^ are not strongly dependent on the initial red- 
shift Zini. To verify that, compare the models Al and A4; 
A2 and A5; A3 and A6. 

However, the value of Zini has influence over the low- 
frequency part of the spectra as can be seen from Figure 7. 
This part of the spectrum is dominated by the population 
of black holes formed at redshifts z > 7. 

It is worth stressing that|de Araujo, Miranda fc Aguiar 



(|2004[) assu ming a Springel and Hernquist ( Springel & Hern 
quist|2003 1 model of star formation obtained a similar result 



for flew- Their spectrum peaks at flaw h « 5 x 10 " at 
i^obs « 200 Hz for a Salpeter IMF. Using h = 0.73 we find 
ficw ~ 9 X 10~^ for their fiducial model. 

This is a factor ~ 5 lower than the maximum ampli- 
tude of flaw obtained by our model A2. However, note that 
p* obtained from 'model 3.0 Gyr' in Figure 2, which corre- 
sponds to model A2 in Table 1, is smaller than the Springel 
and Hernquist CSFR only in the range 4.5 ^ z < 8.2. 
Thus, except for this interval in redshift, the rate of core 
collapse obtained from Springel and Hernquist CSFR is ac- 
tually smaller than that obtained from model A2. 

The cusp in the curves shown in the Figure 7 is pro- 
duced by our choice to the energy flux (see Equations |29| 
and 301. See that the closure energy density (ilcw) is di- 
rectly proportional to the energy flux, and therefore more 
sensitive to its frequency dependence. Here, the specific en- 



ergy fiux is obtained from Equation (331, which takes into 
account the most relevant quasi-normal modes of a rotating 
black hole. 

In particular, we refer the reader to |de Araujo, Miranda| 
& Aguiar ( 2000 1 who discuss the formulation presented here 



and compare it to that used by |Ferrari, Matarrese fc Schnei-| 
|der| ( |1999[ ) where the energy fiux is a function of frequency. 
Thus, their closure energy density is broader than we use 
here. As a consequence, fiow obtained by |Ferrari, Matarrese| 
|fc Schneider] ( |1999[) has a smoother peak than ou rs. However, 
as discussed in de Araujo, Miranda fc Aguiar (20001, both 



formulations presented similar results. 

Since some authors use, instead of Onw. 

1/2 ' 

tional strain S^^ , defined by 



the gr avita- 



Allen fc Romano 



( 1999 1 as 



5h 



3Hi 



47r2 Uoh: 



GW, 



(38) 



we show this quantity in Figure 8. 

A key parameter to determine the values presented in 
Tables 1, 2, and 3 is the efficiency of generation of GWs. 
We take the maximum efficiency found by |Stark fc Piran| 



Table 3. The main results of the models with Afmin = IO^^Mq. 



Model 


^ini 


X (IMF) 


Ts Gyr 




(S/N) 


Cl 


20 


1.35 


1.0 


3.2 


5.4 


C2 


40 


1.35 


1.0 


3.2 


5.6 



^ The model A2 is that which has the smallest values for Hgw- 
As a consequence, from all models presented in Tables 1 — 3, A2 
is that which present the smallest signal-to-noise ratio. 

Concerning for the results presented in Figures 7 and 8 we are 
defining the high-frequency part of the spectra as that for which 
i^obs > 200 Hz. 



stochastic background of gravitational waves 9 




obs 



Figure 7. Spectrum of the gravitational energy density parame- 
ter Oqw Results are shown for the models Al — A6 of the Table 
1. 
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Figure 8. Gravitational strain in Hz Results are shown for 

the models Al - A6 of the Table 1. 



(19861, namely, ecwn 



= 7 X 10~ for an axisymmetric 
collapse resulting in a black hole. 

O n the other hand, more recently, |Fryer, Woosley 

Heger| ( |200l| obtained the efficiency of 2 x 10"^ for a 100 Mq 
black hole remnant. Note that since Claw oc e, if the effi- 
ciency is actually closer to 2 x 10~^, the observed energy 
density in gravitational waves may be divided by a factor of 
35. In this case, of all models here studied only model A6 
will produce (S/N) > 3. 

However, the distribution of ecw in function of the mass 
of a black hole is unknown. In particular, let us think of 
what occurs with other compact objects — namely, the neu- 
tron stars — to see if we can learn something from them. A 
newly born neutron star could lose angular momentum due 
to gravitational waves associated with non-radial oscillations 
I Ferrari, Miniutti fc Pons|2003 1 . This could explain why all 
known young neutron stars are relatively slow rotators. 

The black holes could have had a similar history, i.e, 
they could have been formed rapidly rotating and lost mo- 
mentum to gravitational radiation via their quasi-normal 



modes. If this was the case, the value of ecw could be near 
the maximum one, or in the worst case, it could have a value 
to produce (S/N) > 3 for a LIGO III pair. 



5 SUMMARY AND DISCUSSION 

In this work, we have used the hierarchical formation sce- 
nario derived from the Press-Schechter formalism to build 
the cosmic star formation rate - CSFR in a self-consistent 
way. Our paper differs from earlier works basically in the 
form as is obtained the function (or CSFR). 

In particular, from the hierarchical scenario we obtain 
the baryon accretion rate, ah{t), that supplies the gaseous 
reservoir in the halos. Thus, the term ab(t) is treated as an 
infall term in our model. 

This scenario is in agreement with the cold dark matter 
model of cosmological structure formation, where the first 
sources of light are expected to form in ~ IO^Mq dark mat- 
ter potential at z 20. 

Using we calculate the stochastic background of grav- 
itational waves produced by pre-galactic black holes. We 
show that a significant amount of GWs is produced related 
to the history of CSFR studied here, and this SBGWs can in 
principle be detected by a pair of LIGO III interferometers. 

Note that signal-to-noise ratios (S/N) ~ 90 could be 
obtained if the efficiency of generation of GWs is close to the 
maximum value (eGw,„ax ~ 7 x 10~*), if the IMF produces 
a high number of massive remnants (x = 0.35), and if Zini ~ 
20. Considering a Salpeter IMF (x = 1.35), we obtain signal- 
to-noise ratios (S/N) ~ 10. 

The critical parameter to be constrained in the case of 
a non-detection is ecw. A non-detection would mean that 
the efficiency of GWs during the formation of black holes is 
not high enough. In reality, eGW„ax should be divided by a 
factor > 35 in the case of a non-detection. 

It is worth mentioning that an IMF with x — 2.35 could 
also be responsible for a non-detection same with ecw ~ 
ecw^ax- However, x = 2.35 produces a high number of low 
mass stars that is not in agreement with recent numerical 
simulations of the collapse and fragmentation of primordial 
clouds (see, e.g., Abel, Bryan fc Norman|[2002 1 . 

Another possibility for a non-detection is that the pre- 
galactic stars are such that the black holes formed had 
masses > 5QOM0. In this case, the GW frequency band 
would be out of the LIGO bandwidth. 

However, considering black holes formed from stars with 
masses 25 M0 < m < 140 M0, then the sensitivity of the 
future third generation of detectors could be high enough 
to increase one order of magnitude in the expected value of 
(S/N). Examples of such detectors are the Large Scale Cryo- 
genic Gravitational Wave Telescope (LCGT) and the Euro- 
pean antenna EGO (see Regimbau & de Freitas Pacheco 
1 200 6 and the references therein for a short discussion on 
this subject). 

Specifically, around 650 Hz the planned strain noise for 
EGO will be a factor of ~ 4 higher than that provided for 
advanced LIGO configuration. This could represent a gain 
of a factor ~ 5 — 20 for the value of (S/N) considering two 
interferometers located at the same place (see Regimbau fc 
|de Freitas Pacheco|2006] ). Thus, some models in Tables 1 — 3 
could survive with (S/N) > 3 same with ecw ~ 2 x 10~^. 
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In particular, the detection of a background with sig- 
nificant (S/N) would permit us to obtain the curve (or 
f^Gw) versus Vohs- From it, one can constrain at high red- 
shifts and the gravitational wave efficiency (ecw)- Thus, the 
detection and characterization of a SBGWs could be used 
as a tool for study of the star formation at high redshifts. 

It is worth stressing that several astrophysical sources 
can contribute to the background of gravitational waves, as 
mentioned in the Introduction. In principle, it should be 
possible to distinguish different sources from the detected 
gravitational wave spectrum. That is, from the caractheris- 
tics of the observed curve J^Gw versus t^obs- 

For example, in the present work we have shown that 
cosmological stellar black holes (3 < Mbh/M© < 65), 
formed at Zi^i < 20 — 40, produce a stochastic background 
in the frequency range ~ 10 Hz — 5 kHz. In particular, the 
gravitational wave spectra peak at i/obs ~ 200 Hz. If the 
black hole Population forms at low redshifts (for example, 
Zini < 10), both the frequency where f^Gw peaks and the 
minimum frequency of the spectra will be shifted to greater 
frequencies than those presented here. 

However, the shape of f^Gw does not considerably 
change if we consider the same gravitational wave energy 
power spectrum for the sources. On the other hand, more 
massive stars (m > 200 M©) will shifted the peak of the 
spectra for low frequencies. See for a moment the results of 



Marassi, Schneider & Ferrari ( 2009 1 for black hole remnants 
of Population HI stars with masses 100 — 500 M© . Their 
spectrum peaks at Uohs ~ 2.74 Hz (f^Gw ~ 5 x 10"^"") and 
the maximum frequency of the background is ~ 600 Hz. 

Another example can be seen from the work of |Buo-| 
(20051. The authors studied the gravitational 



nanno et al. 



^'^) while the maxi- 
3 kHz and the spec- 



wave background from all cosmic supernovae. Their fiducial 
model peaks at Uohs = 6Hz {flaw ~ 10~ 
mum frequency of the background is 
trum can extent to very low frequencies {vohe 10"** Hz). 
Thus, in principle, it would be possible to identify the sig- 
natures of different backgrounds if we have the curve How 
versus Vohs over a large range in frequency. 

Last but not least, we refer the reader to the work of 



Kauffmann & Haehnelt ( 2000 1 who present a unified model 



for the evolution of galaxies and quasars. Specifically, these 
authors discuss that gas cooling is not efficient in too massive 
structures and so haloes with circular velocity greater than 
600kms~^ could not form stars. If we take into account their 
results then the upper limit, Mmax, in Equation (10 1 should 
be changed for > 

- 3 with 

ll3^ 



1O"M0. 

We checked all the models described in Tables 1 - 
this new upper limit (Mmax ~ W^^Mq). We verify that the 
amplitude of the CSFR decreases slightly at 2; < 3.5 when 
compared with the results obtained using 
(at z > 3.5 we do not observe any modification in the be- 
haviour of p*). For the models with Mmin = IO'^Mq (IO^^Mq) 
there is only a subtle modification in the final results. In 
particular, the signal-to-noise ratios are ~ 3.9% (~ 4.6%) 
lower than those presented in Table 1 (2). For the models 
with Mniin ~ W^'^Mq wc uote a modification ~ 8.6% in the 
results of the Table 3. However, all models presented in Ta- 
bles 1 — 3 satisfy the "two criteria" , as discussed in Section 
4. That is, same using Mmax = lO^^M© the models produce 
(S/N) > 3 and Xr < 1. 



6 ACKNOWLEDGMENTS 

ESP and ODM thank Claudia V. Rodrigues, Ronaldo E. de 
Souza, Jose C.N. de Araujo, and Jose A. de Freitas Pacheco 
for helpful discussions. The authors would like to thank the 
referee for helpful comments that we feel considerably im- 
proved the paper. ESP was financially supported by the 
Brazilian Agency CAPES, and ODM is partially supported 
by CNPq (grant 305456/2006-7). 



REFERENCES 

Abel T., Bryan G.L., Norman M.L., 2002, Sci, 295, 93 
Allen, B. 1997, in J.-A. Marck, & J.-P. Lasota, eds, Rel- 
ativistic Gravitation and Gravitational Radiation, Cam- 
bridge University Press, Princeton, NJ p. 373 
Allen B., Romano J.D., 1999, Phys. Rev. D, 59, 2001 
Bardeen J.M., Bond J.R., Kaiser N., Szalay A.S., 1986, 
ApJ, 304, 15 

Belczynski K., Kalogera V., Bulik T., 2002, ApJ, 572, 407 
Buonanno A., Sigl G., Raffefi G.G., Janka H.T., Miiller E., 

2005, Phys. Rev. D, 72, 084001 
Carr B.J., 1980, A&A, 89, 6 

Carrol S.M., Press W.H., Turner E.L, 1992, ARA&A, 30, 
499 

Chiosi C, Maeder A., 1986, ARA&A, 24, 329 
Christensen N., 1992, Phys. Rev. D, 46, 5250 
Copi C.J., 1997, Apj, 487, 704 

Daigne F., Olive K.A., Silk J., Stoehr F., Vangioni E., 2006, 
Apj, 647, 773 

de Araujo J.C.N., Miranda O.D., Aguiar O.D., 2000, Phys. 

Rev. D, 61, 12, 124015 
de Araujo J.C.N. , Miranda O.D., Aguiar O.D., 2002, MN- 

RAS, 330, 651 

de Araujo J.C.N., Miranda O.D., Aguiar O.D., 2004, MN- 

RAS, 348, 1373 
de Araujo J.C.N., Miranda O.D., 2005, Phys. Rev. D, 71, 

12, 12703 

de Freitas Pacheco J. A., 1997, Astrop. Physics, 8, 21 
Efstathiou G., Bond J.R., White S.D.M., 1992, MNRAS, 
258, Ip 

Ellison S., Songaila A., Schaye J., Petinni M., 2000, AJ, 
120, 1175 

Ferrari V., Matarrese S., Schneider R., 1999, MNRAS, 303, 
247 

Ferrari V., Miniutti G., Pons J.A., 2003, MNRAS, 342, 629 
Flanagan E.E., 1993, Phys. Rev. D, 48, 2389 
Fryer C.L., Woosley S.E., Heger A., 2001, Apj, 550, 372 
Giovannini M., 2009, preprint (astro-ph:0901.3026) 
Gunn J.E., Peterson B.A., 1965, Apj, 142, 1633 
Heger A., Woosley S.E., 2002, Apj, 567, 532 
Hopkins A.M., 2004, Apj, 615, 209 
Hopkins A.M., 2007, Apj, 654, 1175 
Hulse R.A., Taylor J.H., 1974, ApJ, 191, L59 
Hulse R.A., Taylor J.H., 1975, ApJ, 195, L51 
Hulse R.A., Taylor J.H., 1975, ApJ, 201, L55 
Jenkins A., Frenk C.S., White S.D.M., Colberg J.M., Cole 
S., Evrard A. E., Couchman H., Yoshida N., 2001, MN- 
RAS, 321, 372 
Kauffmann G., Haehnelt M., 2000, MNRAS, 311, 576 
Kroupa P., 2007, preprint (astro-ph:0703124) 



Lukic Z., Heitmann K., Habib S., Bashinsky S., Ricker P. 

M., 2007, ApJ, 671, 1160 
Maggiore M., 2000, Phys. Rep., 331, 283 
Marassi S., Schneider R., Ferrari V., 2009, MNRAS, 398, 

293 

Miranda O.D., de Araujo J.C.N., Aguiar O.D., 2004, Class. 
Quantum Grav., 21, S557 

Peacock J. A., 1999, Cosmological Physics, Cambridge Uni- 
versity Press, Cambrige, 682 

Press W.H., Schcchtcr P., 1974, Apj, 193, 425 

Regimbau T., de Freitas Pacheco J.A., 2006, ApJ, 642, 455 

Salpeter E.E., 1955, Apj, 121, 161 

Salvadori S., Schneider R., Ferrara A., 2007, MNRAS, 381, 

647 

Sandick P., Olive K.A., Daigne F., Vangioni E., 2006, Phys. 
Rev. D, 73, 104024 

Scalo J., Fund. Cosm. Phys., 1986, 11, 1 

Schmidt M., 1959, Apj, 129, 243 

Schmidt M., 1963, Apj, 137, 758 

Sheth R.K., Tormen G., 1999, MNRAS, 308, 119 

Songaila A., Cowie L.L., 1996, AJ, 112, 335 

Springel v., Hernquist L., 2003, MNRAS, 339, 312 

Stark R.F., Piran T., 1986, in R. Ruffini, cd., Proc. Fourth 
Marcel Grossmann Meeting on General Relativity. Else- 
vier Science, Amsterdam, p. 327 

Suwa Y., Takiwaki T., Kotake K., Sato K., 2007, ApJ, 665, 
L43 

Suwa Y., Takiwaki T., Kotake K., Sato K., 2007, PASJ, 59, 
771 

Thorne K. P., 1987, in S.W. Hawking and W. Israel, eds, 
Three Hundred Years of Gravitation, Cambridge Univer- 
sity Press, Cambridge, p. 330 

Venkatesan A., 2000, ApJ, 537, 55 

Wilkins S.M., Trentham N., Hopkins A.M., 2008, MNRAS, 
385, 687 



stochastic background of gravitational waves 



